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ABSTRACT 

Wc present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg 2 
of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of - 10~ 15 erg cm~ 2 s" 1 
(~ 5 x 10~ 15 erg cm -2 s _1 ). In order to identify the correct optical counterparts and to assess the sta- 
tistical significance of the X-ray to optical associations we have used the "likelihood ratio technique" . 
Here we present the identification method and its application to the CFHT I-band and photometric 
catalogs. We were able to associate a candidate optical counterpart to ^90% (626) of the X-ray 
sources, while for the remaining ~10% of the sources we were not able to provide a unique optical 
association due to the faintness of the possible optical counterparts (Lab >25) or to the presence of 
multiple optical sources, with similar likelihoods of being the correct identification, within the XMM- 
Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru 
multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic 
data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis of 
this statistically meaningful sample of X-ray sources reveals that for ~80% of the counterparts there 
is a very good agreement between the spectroscopic classification, the morphological parameters as 
derived from ACS data, and the optical to near infrared colors: the large majority of spectroscopi- 
cally identified broad line AGN (BL AGN) have a point-like morphology on ACS data, blue optical 
colors in color-color diagrams, and an X-ray to optical flux ratio typical of optically selected quasars. 
Conversely, sources classified as narrow line AGN or normal galaxies are on average associated with 
extended optical sources, have significantly redder optical to near infrared colors and span a larger 
range of X-ray to optical flux ratios. However, about 20% of the sources show an apparent mismatch 
between the morphological and spectroscopic classifications. All the "extended" BL AGN lie at red- 
shift <1.5, while the redshift distribution of the full BL AGN population peaks at z^l.5. The most 
likely explanation is that in these objects the nuclear emission is not dominant with respect to the host 
galaxy emission in the observed ACS band. Our analysis also suggests that the Type 2/Type 1 ratio 
decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral 
analysis and the most recent modeling of the X-ray luminosity function evolution. 
Subject headings: surveys — galaxies: active — X-rays: galaxies — X-rays: general — X-rays: diffuse 
background 
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1. INTRODUCTION 

The primar y goals of the "Cosm ic Evolution Survey" 
(COSMOS. IScoville et all l2007al ) is to trace star for- 
mation and nuclear activity along with the mass as- 
sembly history of galaxies as a function of redshift and 
environment. Al though there wer e early theoretical 
suggestions (e.g. iSilk fe Reesl 119981) . it was the tight 
relation observed in local galaxies between the black 
holes mass a nd the velocity dispersion (the Mbh-t rela- 
tion see e.g. iFerrarese fe M crritt 2000l; iGebhardt et al.1 
l2000f ) and the fact that the locally inferred black hole 
mass density appears to be broadly consistent with 
the est imates of the mass accr e ted during the quasar 
phase (IFabian fc Iwasawal 119991; lYu fc Tremaind 2002]: 
Elvis. Risaliti. fc Zamorani 2002; Marc oni et all 120041 ; 



Merloni 2004) that made it clear that BH-driven nuclear 



activity and the assembly of bulge masses are closely 
linked. This realization has led to a large number of 
theoretical studies, suggesting feedback mechanisms to 
explain this fundamental link between the assembly of 
black holes (BH) and the formation of spheroids i n galaxy 
halos (|Di Matteo et al.ll2005l : iMenci et al1l2004l and ref- 
erence therein). 

In this framework, the hard X-ray band is by far the 
cleanest one for studying the history of accretion onto 
black holes in the Universe, being the only band in which 
emission from accretion processes clearly dominates the 
cosmic background. The detailed study of the nature of 
the hard X-ray source population is being pursued by 
complementing deep pencil bea m observations with shal- 
lower, larger area surveys (see iBrandt fc Hasingeij 120051 
for a recent review). Indeed, in the recent years, large 
efforts have been dedicated to the optical to radio char- 
acterization of X-ray sources selected at different X-ray 
depths (see, among the most recent: the XBootes survey 
([Brand et al. Il2006l ) , the Serendi pitous Extragalact ic X- 
ray Source Identification SEX SI feckart et alfeool). th e 
Extended Groth strip EGS (iGeorgakakis et all |200H . 
the HELLAS2XMM survey (|Cocchia et all l2007f )). In 
particular, sizable samples of objects detected at the 
bright X-ray fluxes (> 10~ 14 erg cm~ 2 s _1 ) over an 
area of the order of a few square degrees are needed to 
cover the high luminosity part of the Hubble diagram 
and to obtain, together with samples from narrower and 
deeper pencil-beam observations, a well constrained 
luminosity function with a similar number of sources per 
luminosity decade and redshift bin . The results from 
both deep ([Alexander et all 120031: IBarger et all 120031: 
Giacconi et al. 1 120021 : ISzokolv et al.ll2004t iHasineer et al.l 



2001h and shallow survey s (|Fiore et al.ll2003l : lGreen et al.l 
20041 : ISteffen et al.l l2004f ) have unambiguously unveiled 



a differential evolu tion for the low- and high-luminosity 
AGN population (lUeda et all 120031: IBarger et all | 2005l : 



Silverman et al.|2005t lHasinger. Mivaii fc Schmidtll2005l 

La Franca et a l. 2005). However, in these studies 
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evolution of the high-luminosity tail of the obscured 
AGN luminosity function remains a key parameter still 
to be determined. The best strategy to address this issue 
is to increase the area covered in the hard X-ray band, 
and the corresponding optical-NIR photometric and 
spectroscopic follow-up, down to F2-iofc e v ~ 5 x 10 -15 
erg cm" 2 s _1 , where the bulk of the XRB is produced. 
Another important issue in AGN studies is the deter- 
mination of the global Spectral Energy Distributions 
(SEDs), over the widest possible frequency range, of 
different types of AGN. Large samples of AGN have 
been assembled using different selection criteria at 
different frequencies (e.g. radio, infrared, optical-UV, 
X-ray); however, the lack of complete multi- wavelength 
coverage for many of these samples makes it difficult 
the comparison of the properties of different classes of 
AGN. For example, while the average SED of optically 
and radio selected bright quasars ha s been reasonabl y 
well known for more than a decade (|Elvis 
little is still known about the SEDs of lower luminosity 
X-ray selected AGNs and, in particular, of the obscured 
ones. This reflects into a significant uncertainty in 
the estimate of the bolomctric correction (kboi, i.e. 
the correction from the X-ray to the bolomctric lu- 
minosity) to be applied to the observed luminosity of 
X-ray selected AGN and, in turn, on the estimates 
of the Eddington lumi nosities and of the masses of 
the central black holes (jFabianl I2004T ) . Recent results 
on high-redshift quasars have shown that the overall 
X-ray to optical spectral slope (a x, usually measured 
between 2500 A and 2 keV, e.g. iTananbaum et al.lll979l) 
and the corresponding fc& ;, is a function of the AGN 
luminosity, with the low luminosity objects having a 
lower khni (IVignah. Brandt fc Schneiderl 12001 IFabian! 
12004 ISteffen et al.ll2006l ). It is therefore clear that, in 
order to get a complete census of the luminosity output 
of AGN, a multi- wavelength project is needed: only in 
this way can the SEDs from radio to X-rays of a large 
sample of AGN, selected at different frequencies, be 
compiled. 

The XMM-Newton wide-field survey in the COSMOS 
field (hereinafter: XMM-COSMOS) is an important 
step forward in addressing the topics described above. 
The ~ 2 deg 2 a rea of the HST/ ACS COSMOS Trea- 
sury program (|Scoville et all l2007bD . bounded by 
9 h 57.5 m <R.A.< 10 h 03.5 m ; l d 27.5 m <DEC< 2 d 57.5 m , 
has been surveyed with XMM-Newton for a total of 
~800 ks during A03, and additional 600 ks have been 
already granted in A04 to this project. When com- 
pleted, XMM-COSMOS will provide an unprecedently 
large sample of X-ray sources (> 2000), detected on 
a large, contiguous area, with complete ultraviolet 
to mid-infrared (including Spitzer data) and radio 
coverage, and almost complete s pectroscopic follq w-up 
granted through t he zCOSMOS (iLillv et all l2007f) and 
Magellan/IMACS (jlmpev et al.ll2007f ) projects. 
The XMM-COSMO S project is described in 
lHasinger et alj (|2007| ). The X-ray point source counts 
from the first 800 ks of the X MM-Newton observation s 
obtained so far are presented in lCappelluti et all (|2007h . 
while the pr operties of extend e d and diffuse sources are 
described in lFinoguenov et al.l (|2007). First results from 
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the spectral analysis of AGN with know n redshift and 
high counting statistics are presented in Mainieri et al. 

pool . 

This paper presents the optical identification of X- 
ray point sources detected in this first year of the 
XMM-Newton data, and the first results on the mul- 
tiwavclcngth properties of this large sample of X-ray 
selected AGN. 

The paper is organized as follows: Section 2 presents 
the multiwavelcngth datasets and describes the method 
used to identify the X-ray sources and its statistical 
reliability; the X-ray to optical and near infrared prop- 
erties are presented in Sect. 3, along with preliminary 
results from the ongoing spectroscopic follow-up; in 
Sections 4 and 5 we discuss and summarize the most 
important results. Throughout the paper, we adopt 
the cosmological para meters Hp = 70 km s _1 Mpc -1 , 
O m =0.3 and ft A =0.7 (jSpergel et al.ll2003fh 

2. X-RAY SOURCE IDENTIFICATION 

2.1. Optical and X-ray datasets 

The XMM-COSMOS X-ray source catalog comprises 
1390 different point-like X-ray sources detected over 
an area of ^2 deg 2 in the first-ye ar XMM-Newton ob- 
servations of the COS MOS field (jHasinger et all 120071 ; 
ICappelluti et all l2007[ K In this paper, we limit our 
analysis to the sources detected in the first 12 XMM- 
COSMOS observations (over a total of ~1.3 deg 2 ), for 
which both the X-ray source list and the optical catalogs 
were in place in June 2005, and a substantial spectro- 
scopic follow -up already exists. Th ese fields are flagged 
in Table 1 in lHasinger et all (|2007f ). and the average ex- 
posure ti me is of ~ 25fcs (~ 4 ks in the overlapping re- 
gion, scc lHasinger et aLll2007l for more details). The cat- 
alog paper with the identification from the entire XMM- 
COSMOS survey is in preparation and will rely strongly 
on the identification procedure described here. 
The 12 fields X-ray catalog (hereinafter: 12F catalog) 
comprises 715 X-ray sources detected over a total of 
~1.3 deg 2 . From the 12F catalog we removed 20 sources 
classified as extended by the detection algorithm (see 
ICappelluti et alj 120071 for a description of the adopted 
detection algorithm). The observed X-ray emission 
from most of these sources is likely to be due either to 
groups/clusters of galaxies, or to the contribution of two 
or more X-ray sources close to each other. In both cases, 
an association with a unique optical counterpart is not 
possible. The total number of point-like X-ray sources 
in the 12F catalog, detected in at least one of the X-ray 
bands is therefore 695. Of these, 656 are detected in the 
soft band (0.5-2 keV), 312 in the medium band (2-4.5 
keV; 38 only in the medium band), and 47 in the hard 
band (4.5-10 keV; 1 only in the hard). The X-ray cen- 
troids of these 695 point-like X-ray sources have been as- 
trometricall y calibrated using t he SA S task eposcorr, as 
described in ICappelluti et all (|2007f ): the resulting shift 
of ~1" (A(a)=0.99"; A(S)=0") was applied to all source 
positions. 

As a first step in the identificati on process, we used the I - 
band CFHT/Megacam catalog (|Mc Cracken et al.ll2007f h 
which, although slightly shallower than other available 
dat a (e.g. the Subaru B,g,V,R,I and z photometric data, 
see lCapak et al.|[2007l ). has the advantage of having reli- 



able photometry even at bright magnitudes (Iab ?S 19), 
where the Subaru photometry starts to be significantly 
affected by uncertainties due to saturation. 
At magnitudes brighter than Iab = 16, only ~0.24 
sources are expected by chance in a 3" error-circle around 
all the 695 X-ray sources on the basis of the background 
counts of the Megacam catalog; therefore we considered 
as secure optical identifications all the 21 sources brighter 
than I^s=16 in the CFHT catalog and within 3" from 
the X-ray centroids. For all the remaining X-ray sources 
(674) we used the method described in Sect. 12.21 
Then, we have cross-correlated the optical counterparts 
with t he multicolor phot o- z catalog (June 22th 2005 re- 
lease, iCaoak et al.ll2007l) 2 . All the magnitudes are in 
the AB system Ij0kdll971l) . if not otherwise stated. 



2.2. The method 

Typical error-circles of XMM-JVeroion data (~ 5" arcscc 
radius at 95% confidence level 22 ) often contain more than 
one source in deep optical images, so that the identifica- 
tion process is not always straightforward. Obviously, 
this is particularly true when the candidate counterparts 
are faint. We therefore decided to use the "likelihood 
ratio" {LR) technique, i n order to properly identify 
the optical counterparts ([Sutherland fc Saunders! 1 19921 : 
ICiliegi et al.ll2003l : iBrusa et al.l l2005fh The LR is defined 
as the ratio between the probability that the source 
is the correct identification and the corresponding 
probability of being a back ground, unrelated object 
([Sutherland fc Saunders! 1 1992f ). i.e.: 
LR = « m )f<r) 

n{m) 

where q(m) is the expected probability distribution, 
as a function of magnitude, of the true counterparts, 
f(r) is the probability distribution function of the 
positional errors of the X-ray sources assumed to be 
a two-dimensional Gaussian, and n(m) is the surface 
density of background objects with magnitude m. For 
the calculation of the LR paramet ers we followed an d 
improved the procedure described in lBrusa et al.l (|2005f ). 
For the f(r) calculation, we used the statistical error 
as computed from the detection procedure and tested 
against the pattern of observati ons with extensive Mont e 
Carlo simulations presented by ICappelluti et"aTI (|2007f) . 
We also added in quadrature a systematic error of 0.75": 
as noticed by Loaring et al. (2005), this additional 
component may be due to residual uncertainties in the 
detector geometry, and may represent a fundamental 
limit to the positional accuracy of XMM-Newton. We 
adopted a 3" radius for the estimate of q(m), obtained 
by subtracting the expected number of background 
objects {n(m)) from the observed total number of 
objects listed in the catalog around the positions of the 
X-ray sources. Since on the basis of several results in the 

20 The new version of the photometric catalog has been released 
at the time of the submission of this paper. The (small) differences 
between the two photometric catalogs do not affect the results 
presented in this work. 

21 We adopted the following AB-Vega conversion: 
I(Vega)=I(AB)-0.4, R-K(Vega)=R-K(AB)+1.65. 

22 Such a radius represents the radius for which ~95% 
of the XMM-Newton sources in the SSC catalog are 
associated with USNO A. 2 sources (see Fig. 7.5 in 
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Fig. 1. — Observed magnitude distribution of all the optical objects detected in the I band within a radius of 3" around each of the 674 
point-like X-ray sources fainter than Iab=16 (solid histogram), together with the expected distribution of background objects in the same 
area (n(m), dashed histogram). The difference between these two distributions (dotted histogram) is the expected magnitude distribution 
of the optical counterparts. The smooth curve fitted to this histogram (dotted line) has been used as input in the likelihood calculation 
(q(m)). 



literature (see e.g. iFiore et al.1 120031 : iDella Ceca et al.l 
[200l lLoaring et all 120051) . a large fraction of the possi- 
ble counterparts are expected to be found within a 3" 
radius, this choice maximizes the statistical significance 
of the over-density around the X-ray centroids, due to 
the presence of the optical counterparts. With this 
procedure, q(m) is well defined up to Iab ~23.0-23.5. At 
fainter magnitudes, the number of objects in the error 
boxes around the X-ray sources turns out to be smaller 
than that expected from the field global counts n(m). 
Formally, this would produce an unphysical negative 
q(m), which, in turn, would not allow the application 
of this procedure at these magnitudes. The reason for 
this effect is the presence of a large number of relatively 
bright optical counterparts (Iab = 16-23) close to the 
X-ray centroids. These objects, which occupy a non- 
negligible fraction of the total area of the X-ray error 
boxes, make it difficult to detect fainter background 
objects in the same area. As a consequence, the n(m) 
estimated from the global field is an overestimate of 
the observed n(m) at faint magnitudes in the X-ray 
error boxes. In order to estimate the correct n(m) to 
be used at faint magnitudes in the likelihood calcula- 
tion, we have randomly extracted from the Megacam 
catalog 1000 optical sources with the same expected 



magnitude distribution of the X-ray sources, and we 
computed the background surface density around these 
objects. As expected, we found that, indeed, the n(m) 
computed in this way is consistent with the global n(m) 
at Iab<23.0, but is significantly smaller than it (and 
smaller than the observed counts in the error boxes) at 
fainter magnitudes. Therefore, the input n(m) in the 
likelihood procedure was the global one for Iab < 23 and 
that derived with this analysis around sources with the 
same magnitude distribution as the "bright" optical 
counterparts for Iab >23.0. This allowed us to identify 
a few tens of very faint sources that would have been 
missed without this correction in the expected n(m). 

Fig. 1 shows the observed magnitude distribution of all 
optical objects present in the / band catalog within a 
radius of 3" around each X-ray source (solid histogram) , 
together with the expected distribution of background 
objects in the same area, estimated using the procedure 
described above (dashed histogram). The difference 
between these two distributions (dotted histogram) 
is the expected magnitude distribution of the optical 
counterparts. The smooth curve fitted to this histogram 
(dotted curve) has been used as the input in the 
likelihood calculation q(m); the normalization of this 
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TABLE 1 

Summary of optical identification of point-like X-ray sources 





all 


identified 


unidentified 


Sources with I< 16 (I bright) 


21 


21 




Sources with I> 16 


674 


587 


87 


Additional sources identified with K 




46 




Total 


695 


654 


41 






rcliable a ambiguous 15 








626 28 





a We classified as "reliable" all the 21 I-bright sources, the 46 sources identified from 
the K-band, and 559 objects for which there is only one object with LR> L-Rth or 
the ratio between the highest and the second highest LR value in the I and/or in the 
K-band is greater than 3. See Sect. 2.3 for details. b We classified as "ambiguous" 
all the sources for which the ratio of LR values of the possible counterparts in the 
I-band and in the K-band is smaller than 3. See Sect. 2.3 for details. 
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Fig. 2. — (Upper panel): examples of I-band finding charts 
of sources with a unique optical counterpart or securely identified 
from the likelihood analysis. (Middle panel): examples of I-band 
finding charts of sources with ambiguous counterparts. (Lower 
panel): examples of K— band finding charts of sources identified 
from the K-band catalog. In all the panels, the circle is centered 
on the XMM— Newton position and has a radius of 5". 

curve has been set to 0.76, corresponding to the ratio 
between the integral of the q(m) distribution and the 
total number of the X-ray sources. For the threshold 
value in the likelihood ratio we adopted Li? t h=0.4, 
which maximizes the sum of sa mple reliability and com- 
pleteness in the present sample ([Sutherland fc Saundersl 
fl99l lBrusa"et~al1l2005h . 



2.3. Results 

We computed the likelihood ratio for all the optical 
sources within 7" of the 674 X-ray centroids (a total of 
2158 sources), finding 730 optical sources with LR>LR t h- 
The expected number of true identifications among these 
objects, computed by summing the reliability of each 
of them, is ^524 (~78% of the total number of X-ray 
sources). The number of X-ray sources which have at 
least one optical source with LR>LR t h is 587 (87% of 
the sample). The remaining 87 X-ray sources either are 
associated with LR<LR t h counterparts (74) or do not 



have any optical counterpart in the catalog within 7" 
from the X-ray centroids (13). Wc have therefore visu- 
ally checked these 87 unidentified sources. In 24 cases we 
found a clear, relatively bright optical source within the 
X-ray error-circle in the I-band image. These objects 
were missing in the input CFHT catalog used to com- 
pute the likelihood ratio (mainly because close to bright 
objects and/or in CCD masked regions) or they were 
present but incorrectly associated with a much fainter 
magnitude. The I-band magnitudes of these 24 sources 
were therefore retrieved from the Subaru multicolor cat- 
alog, and added to our list of possible counterparts. It 
is important to note that this problem affects only a mi- 
nority of the X-ray source counterparts (~3.5%, most 
of them close to bright stars and/or defects in the op- 
tical band images) and it can be easily solved with the 
adopted procedure. 

We then ran the likelihood analysis again, this time us- 
ing as input catalog the K-band (KPNO /CTIO) data 
extra cted from the multicolor catalog (see ICapak et alJ 
l2007h . The main advantage of using also this near- 
infrared catalog is due to the fact that the X-ray to near- 
infrared correlation for AGN is much tighter than the one 
in th e optical bands (jMainieri et al.l l2002t iBrusa et al.1 
l2005f ). Moreover, although the currently available Ri- 
band data are shallower than the optical ones, their use 
is potentially important to find the reddest optical coun- 
terparts, which are a not negligible fract ion of the identi- 
ficati ons of faint X-ray sources (see, e.g. JAlexander etHTI 
l200lh . 

Indeed, using this catalog and adopting the same 
value of LR t h=0.4, we found likely counterparts for 46 
of the 87 unidentified sources, that include all the 24 
discussed above and 22 additional red, faint sources 
with unambiguous identification in the K-band (see also 
iMignoli et afll2004h . 

Summarizing, on the basis of the likelihood analysis 
applied to the I-band and K-band catalogs, we have pos- 
sible counterparts for 633 out of 674 X-ray sources (654 
out of 695 when the 21 sources brighter than 1=16 are 
included). We have further tentatively divided the 654 
proposed identifications in "reliable" and "ambiguous" . 
In the first class we include the 21 sources with Iab<16, 
the 559 sources for which either there is only one object 
with LR> LR t h or the ratio between the highest and 
the second highest LR value in the I or K-band is 
greater than 3, and the 46 sources identified only in the 
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Fig. 3. — (Upper panel): Histogram of the distances between 
the X— ray and optical counterparts for all the 654 sources identified 
as discussed in Sect. 2.3 (solid histogram) and the subsample of 
sources with spectroscopic redshift and classification (shaded his- 
togram). The inset in the top-right part of the plot shows the 
cumulative distribution of X-ray to optical distances for all the 
sources (bold curve) and for the subsample of sources with spec- 
troscopic redshift (dashed curve). The dashed line at 3" marks 
the radius which encloses 90% of the X-ray counterparts (see dis- 
cussion in Section 2.2). (Lower panel): Histogram of the I-band 
magnitudes distribution for the 633 primary identifications fainter 
than Iab=16. The solid and shaded histograms have the same 
meaning as in the upper panel. The dashed histogram shows the 
lower limit to the Iab magnitudes, corresponding to the brightest 
objects in the error-circles, for the 41 X-ray sources without an 
optical identification. The inset in the top-right part of the plot 
shows the cumulative distribution of X— ray to optical distances for 
all the sources with 1^,16 (bold curve) and for the subsample of 
sources with spectroscopic redshift (dashed curve). 

K-band. The total number of objects in this class is 
therefore 626. The combined use of optical and NIR 
catalogs for the identification process led to a number 
of "reliable" counterparts larger than that inferred from 
the use of the optical catalog only (524). However, we 
note that ~ 5% of these reliable counterparts can still 



be misidentifications. The other 28 sources comprise all 
the sources for which the ratio of LR values in the I- or 
K-band of the possible counterparts is smaller than 3, 
and they are considered as ambiguous. In the following 
we tentatively identify these 28 X-ray sources with the 
optical counterparts with the highest LR value. 
Tabic 1 gives a summary of the identifications, while 
Figure 2 shows a few examples of finding charts for 
sources representative of the different types discussed 
above. 

For the remaining 41 sources, the possible counterparts 
have on average faint or very faint optical magnitudes, 
none of them has LR> LR t h, and most of them lie 
at A(X — O) >2". For a fraction of them the true 
counterpart of the X-ray source may indeed be very 
faint and/or undetecte d in the optical bands (see e.g. 
iKoekemoer et alll2004i ). Another possibility is that the 
X-ray emission originates from a group of galaxies at 
high-redshift with optically faint members: in this case, 
the X-ray extension of the putative group emission may 
be comparable to or smaller than the XMM-Newton PSF 
and the source can have been classified as point-like. Fi- 
nally, a small fraction of these objects can be spurious 
X-ray sources, as suggested by the fact that the percent- 
age of sources wit h low X-ray detection likelihood (e.g. 
DETML < 9; see ICappelluti et alj|2007l for a definition 
of DETML) is significantly higher for these unidentified 
sources (~ 30%) than for the optically identified X-ray 
sources (~ 10%). A detection with Chandra with signif- 
icantly smaller error-circles, would definitively discrimi- 
nate among these different possibilities. 

2.4. Statistical properties 

The distribution of the X-ray to optical distances and 
of the I-band magnitudes for the X-ray sources identi- 
fied as described in Sect. 2.3 are shown in Fig. 3 (solid 
histograms). For the X-ray sources with more than one 
optical counterpart with LR>LR t h, the X-ray to optical 
distance and the I band magnitude of the optical ob- 
ject with the highest LR is plotted. In the lower panel 
the dashed histogram shows the lower limits to the Iab 
magnitudes, corresponding to the brightest objects in the 
error circles, when the 41 X-ray sources without an op- 
tical identification are also added. 

About 90% of the reliable optical counterparts have an 
X-ray to optical separation (A(X — O)) smaller than 3" 
(see inset in the upper panel of Fig. 3), which is consis- 
tent with or even better than the ~ 80% within 3" found 
in XMM-Newton data of comparable depth (see e.g. 
Fiore et all 120031 : iDella Ceca et all 120041 : lLoaring et all 
2005h . The improvement is likely given by the combina- 
tion of the accurate astrometry of the p ositions of the X- 
ray s ources (tested with simulations in ICappelluti et all 
l2007f ). and the full exploitation of the multiwavelength 
catalog (optical + K) for the identification process (that 
makes "reliable" a larger number of sources). In addi- 
tion, because of the geometry of the final mosaic, most 
of the X-ray sources are closer than 5' to the center of 
at least one pointing, so that the corresponding PSF is 
not significantly deteriorated. The majority of the X-ray 
sources (~73%) have counterparts with Iab magnitudes 
in the range 20 < Tab <24, with two tails at fainter 
(~12%) and brighter (~15%) magnitudes. The median 
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Fig. 4. — Reds hift distribution for th e sample of 243 sources with "secure" spectr oscopic information from Magcllan/IMACS data 
IITrump et al.ll2007ft '), zCOSMOS project ILillv et al.ll2007l) and from the SDSS archive lAdelman-McCarthv eTaI|[2u05l) . The black filled 
histogram shows the contribution of spectroscopically classified Broad Line AGN. The redshift windows in which each broad line (Ha, H/3, 
Mgll, CIII], CIV) is visible in the IMACS/zCOSMOS spectra (wavelength range ~5400-9000 A) are also reported in the upper part of the 
Figure. 



magnitude of optical counterparts is Iab =21.65. 
These distributions do not change significantly if we 
make a different choice for the optical counterpart of the 
28 "ambiguous" sources, i.e. the second most likely op- 
tical counterpart is considered. Therefore, from a sta- 
tistical point of view the properties of these 654 X-ray 
sources (sec Table 1) can be considered representative of 
the overall X-ray point source population, at the sampled 
X-ray fluxes. 

2.5. Redshift distribution 

Spectroscopic redshifts for the proposed counterparts 
are available from the Magellan /IM AC S observation 
camp aign (193 objects, llmpey et al.l 120071 iTrump et al.1 
12007b a nd the first run of zCOSMOS observations (25 
objects. ILillv et al"1l2007t J. P. Kneib, private communi- 



cation) and/or were al ready present in the SPSS s urvey 
catalog (48 objects, lAdelman-McCarthv et all 120051 : 
iKaufmann et al.ll2003r 3 ). for a total of 248 independent 
redshifts (five of them of galactic stars). Spectroscopic 
information and classification is therefore available for a 
substantial fraction of our optical/infrared counterparts 
(248/654, -38%). The X-ray to optical distances and 
the I-band magnitude distributions of the subsample 
with spectroscopic redshift are shown as shaded his- 
tograms in both panels of Fig. 3. The spectroscopic 
completeness is ~50% at Iab <22 and drops to ~35% 
at 22< Iab <23 and to ~18% at 23< I A b <24. 
While a more detailed analysis of the optical spectra 

23 These sources have been retrieved from the NED, NASA Ex- 
tragalactic Database 
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Fig. 5. — Redshift versus I-band magnitude for sources spectro- 
scopically classified as BL AGN (green circles) and NOT BL AGN 
(yellow squares). 
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Fig. 6. — Photometric redshifts versus spectroscopic redshifts 
for the sources with secure redshifts in the 12F sample. Red circles 
mark extended sources, blue triangles mark point-like objects, on 
the basis of the ACS data. The dotted lines limit the locus for 
which \z spec - z phot \ < 0.15 X (1 + z spcc ) 



will be presented in future papers, based also on more 
data that are rapidly accumulating from the on-going 
projects, for the purposes of the present paper we 
divided the extragala ctic sources ( 243 o bjects) in only 
two classes, following iFiore et al.1 (|2003[ ): objects with 
broad optical emission lines (BL AGN, 126 in total) and 
sources without a clear, broad emission signature in the 
optical band, but with narrow emission or absorption 
lines (NOT BL AGN, 117). The latter class is therefore 
a mixed bag of different types of objects, and comprises 
mainly narrow-lined AGN, low-luminosity AGN, star- 
bursts and normal galaxies. Figure 4 shows the redshift 
distribution for all the extragalactic sources in the spec- 
troscopic sample (solid empty histogram) along with 
the redshift distribution of the sources spectroscopically 
classified as BL AGN (filled histogram). At low redshift 
(z<l) the spectroscopic identifications are dominated 
by objects associated with NOT BL AGN, while at 
high redshifts (z>1.5) almost exclusively B L AGN are 
detected. As a lready noted bv lEckart et al.l (|2006[ ) and 
ICocchia et all (|2007l) . the small number of narrow line 
objects at z>1.5 is most likely due, at least in part, 
to selection effects: high-redshift narrow line AGN are 
on average fainter in the I-band than BL AGN and 
a significant fraction of them is expected to be below 
the optical magnitude limit for spectroscopic follow-up. 
This is shown in Fig. 5, where the spectroscopic redshifts 
are plotted versus the I-band magnitudes of the optical 
counterparts, for both BL AGN (green circles) and NOT 
BL AGN (yellow squares): while BL AGN show a large 
spread in I magnitudes over a large redshift range, the 
distribution of NOT BL AGN in this plane is much 
narrower and shows a rapid rise in I magnitu de with 
redshift (see also Fig. 14 in iTrump et al.l 120071 paper) . 
Conversely, the small number of BL AGN at low-z can 
be ascribed to several factors: 1) spectroscopic mis- 
classification: given the wavelength range (^5400-9000 
A) of both the IMACS and the VIMOS zCOSMOS 
data, at low redshift (z <1) broad emission lines are 
not well sampled (e.g. Ha is not visible at z> 0.4, H/3 
is not visible at z> 0.8 and Mgll enters the spectral 
window at z~l), while narrow lines (e.g. [Oil]) and/or 
absorption features (e.g. Balmer break) can be more 
easily detected; 2) especially in the case of moderately 
weak Seyferts, the relative contribution of the host 
galaxy and the AGN can be such that the broad lines 
are diluted in the stella r light and can be revealed only 
in hi gh S/N spectra flMoran. Filippenko fc Chornockl 
[2rMlSe^er"gnini et al.ll2003l) . 

To obtain a complete redshift distribution for the X- 
ray sources, photometric redshifts should be computed 
for the sources without spectroscopic data. As a train- 



ing sa mple, we used the code described in lBender et al.l 
(|2001h and tested it on the sample of objects with se- 
cure spectroscopic redshifts. For objects classified as 
"extended" in the ACS images (see Section 3.1), we ap- 
plied the same semi-empirical templates of non-active 
galaxies that have been used to compute the photomet- 
ric re dshift in the Fornax Deep Fields (see Gabasch et al. 
120041 and reference s therein) and GOODS-South field 
(|Salvato et al.ll2~006h : in both applications, a <t(Az/(1 + 
z)) ~0.05 has been obtained. For objects that have been 




Fig. 7. — Full Width Half Maximum (FWHM) vs. ACS I-band magnitude for all the sources (small black symbols) in the ACS photometric 
catalog. Overplotted as blue triangles (red circles) are the X-ray sources classified as point-like (extended) from ACS. The histograms on 
the right part of the figure show the distribution of FWHM for the sources detected only in the soft band (dashed histogram) and for the 
sources detected only in the hard band (solid histogram). Both histograms are normalized to the total number of sources in each class. 



classified as point-like sources, only AGN te mplates have 
been adopted (see also lMainieri et al"1l2005l and reference 
therein). In Fig. 6 our best attempt so far is shown. 
Red circles are optically extended X-ray selected sources, 
while blue triangles represent point-like objects. A rela- 
tively high fraction of both extended (~73%) and point- 
like (~66%) sources have \z spec — z p hot\ < 0.15 x (1 + 

^spec) • 

3. MULTIWAVELENGTH PROPERTIES 

3.1. ACS morphology: stellar vs. extended objects 

The catalog of the primary optical counterparts of the 
X-ray sources was then cross-co rrelated with the June 
2005 version of the AC S catalog (jLeauthaud et al.ll2007l : 
iKoekemoer et al.ll2007t ). from the first HST cycle (Cy- 
cle 12), in order to gather some preliminary information 
on the morphological classification of the X-ray sources. 
Since the available ACS catalog does not cover the en- 
tire XMM- COSMOS area analyzed here, wc could use 



ACS information only for 524 optical sources, out of the 
total of 654. For these sources, we retrieved from the 
catalog the measure of the Full Width Half Maximum 
(FWHM, in image pixels), the ACS I-band magnitude 
and a parameter that defines the morphological classifi- 
cation (stellar or extended) on the basis of the analysi s 
of the available data. Following iLeauthaud et aTl (|2007f ). 
objects were divided in "stellar/point-likc" (hereinafter: 
point-like) and "optically extended" ' (hereinafter: ex- 
tended) on the basis of their position in the plane de- 
fined by the peak surface brightness above the back- 
ground level and the total magnitude (see Fig. 6 in 
ILeauthaud et al.ll2007t ). The FWHM versus the I-band 
magnitude is shown in Fig. 7; in this figure, the blue tri- 
angles indicate sources classified as "point-like" , and red 
circles sources classified as "extended". More than 50% 
of the primary IDs have stellar (or almost stellar; FWHM 
< 3 pixels) profile on ACS data. This is particularly true 
for the very soft sources (i.e. detected only in the soft 
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Fig. 8. — U-B vs. B-V diagram for all the fields objects with B> 19 (small black poi nts) classified as point-like in the full ACS catalog 
and with an error in all the three optical bands smaller than < 0.15 (Capak et al. 2007). The locus occupied by stars is clearly defined in 
this diagram, with the two densely populated regions in the blue (upper left of the black points envelope; hot subdwarf stars) and red (lower 
right; dwarf M stars) parts of the sequence. The horizontal solid line marks the color cut for the selection of ultraviolet excess objects. 
Overplotted as blue triangles are the counterparts of X— ray sources classified as point-like from ACS, which satisfy the same selection in 
the photometric errors. Spectroscopically confirmed BL AGN and NOT BL AGN are also indicated (green and yellow symbols overplotted 
on the blue ones, respectively). Finally, the expected tracks of quasars from redshift to ~3.5 in this color-color diagram is also reported 
(see text for further details). 



band, see dashed histogram in Fig. 7). The situation is 
completely different for the really hard sources (i.e. no 
detection in the soft band, solid histogram in Fig. 7), for 
which most of the counterparts (~80%) are associated 
with extended sources. 

A subsample of 214 objects (108 BL AGN, 101 NOT 
BL AGN, 5 stars) has both spectroscopic information 
and ACS classification. In agreement with the expecta- 
tions, we find a good correspondence between ACS and 
spectroscopic classification: the majority of BL AGN are 
classified as point-like by ACS (86/108, -80%), while the 
majority of NOT BL AGN are classified as extended by 
ACS (85/101, ~85%). A total of 38 sources show a "mis- 
match" between the morphological and spectroscopic 
classification and we will come back to these sources later 



in the paper. 

3.2. Optical color-color diagrams for point-like objects 

Figure 8 shows the U-B vs. B-V diagram for all field 
objects (small black points) with B>19 classified as 
point-like in the full ACS catalog and with an error in 
all the three optica l bands smaller than 0.15 magnitude 
(|Caoak et alj 12007). The locus occupied by stars is 
very well defined in this diagram, with the two densely 
populated regions in the blue (upper left of the black 
points envelope; hot subdwarf stars) and red (lower 
right; dwarf M stars) parts of the sequence. Most of 
the points at U-B<0.3 (UVX objects) are expected to 
be AGN at z<2.3, while the objects on the right of 
the stellar sequence are likely to be compact galaxies 
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Fig. 9. — The I-band magnitude (Vega system) vs. the soft (0.5-2 keV) flux for all the primary counterparts in the sample (all symbols), 
and the subsamples of objects classified as point-like (blue circles) and extended (red circles) from ACS morphology. The shaded area 
represents the region typically occupied by known AGN (e.g. quasars and Seyferts) along the correlation log(X/0) = 0±1. Spectroscopically 
confirmed BL AGN and NOT BL AGN are also indicated (green and yellow symbols, respectively). Asterisks mark the objects tentatively 
identified with stars (see Sect. 3.2 and 3.3). Lower limits to the I band magnitude correspond to the magnitude of the brightest objects in 
the error-circles of the unidentified X-ray sources. 



classified as point-like. Ovcrplotted as blue triangles 
are the X-ray sources classified as point-like from ACS 
which satisfy the same selection in the photometric 
errors. Spectroscopically confirmed BL AGN and 
NOT BL AGN are also indicated (following Fig. 5, 
green and yellow symbols overplotted on the blue ones, 
respectively). Finally, the expected track of quasars 
from redshift to 3.5 in this color-color diagram is also 
reported, computed us ing the SP SS AGN template f rom 
(jBudavari et all l200lt see also ICapak et all 120071 for 
further details). 

It is quite reassuring that the majority (> 80%) of X-ray 
selected point-like quasars occupy the classical QSO 
locus and would have been selected as outliers from 
the stellar locus in this color-color plot. A few (5 to 8) 
sources are within or close to the locus of "red" stars. 



So far we have spectroscopic data only for two of these 
objects and both of them are stars. All these sources 
are detected only in the soft X-ray band and show a 
low (<0.1) X-ray to optical flux ratio (X/O) 24 and are 
therefore likely to be stars. A similar number (8 to 10) 
of X-ray sources are within or close to the locus of blue 
stars. The spectroscopic identifications available for 
these objects include three BL AGN at z> 1.4; a fourth 
object is likely to be a star on the basis of the bright I 
band magnitude and low X/O ratio. 

3.3. X-ray to optical flux ratio 

24 The I— band flux is computed by converting I magnitudes into 
monochrom atic fluxes and then multiplying them by the width of 
the I filter (Zombcck 1990). 
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Fig. 10. — R-K color vs. K band magnitude (Vega system) for the X-ray sources counterparts with ACS information. Symbols as in 
Fig. 9. Sources spectroscopically identified with BL AGN (green) and NOT BL AGN (yellow) are also highlighted. The histogram on the 
right part of the figure shows the distribution of R-K colors for the sources with HR> —0.3 (i.e. candidate X-ray obscured AGN, solid 
histogram) and the sources detected only in the soft band (HR=— 1, dashed histogram). 



Another independent check on the agreement between 
optical identification, morphological analysis and spec- 
troscopic breakdown is shown in Fig. 9, where the 
I-band magnitude (Vega system) is plotted versus the 
soft X-ray flux for all the X-ray sources (black empty 
symbols), ACS-point-like objects (blue symbols) and 
ACS-extended sources (red symbols). 
As discussed in Sect 2.1, the sources brighter than 
Iab=16 are identified "a priori": all of them have been 
visually inspected and the majority of them turned out 
to be associated with bright stars and have point-like 
morphology from ACS. A sample of 8 additional stars 
has been tentatively identified on the basis of the 
information from U— B vs. B— V diagram, the bright 
optical magnitudes and the low X-ray to optical flux 
ratio. All the sources identified with stars (a total 
of 25) are plotted as asterisks in Fig. 9 and they 
constitute ~4% of the soft X-ray selected sample. This 



value is similar to that repo rted in the ROSAT Deep 
Survey (|Le hmanu et al. 200l]) and in the ChaMP survey 
(| Green et al.l [20040 , from the optical identification of 
large samples (100-300) of soft X-ray selected sources 
at limiting fluxes similar to that of COSMOS. 

If objects associated with stars are removed, point-like 
sources occupy preferentially the locus in this plane with 
X/O in the range 0.1-^10, where most of the broad line 
quasars from previous optical and soft X-ray selected 
surveys lie (hatched region in Fig. 9). Also most of the 
extended sources are within the same range of X/O ra- 
tio, but with more significant tails toward both low and 
high values of X/O. For most of these sources, thanks to 
the superb ACS resolution, it will be possible to resolve 
the nucleus and the host galaxy. Optically bright (i.e. 
X/O<0.1) extended sources are preferentially identified 
with nearby galaxies, in which the optical luminosity is 
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mainly due to the integrated stellar light; for the majority 
of them, the high-level of observed X-ray flux suggests 
that some activity is taking place in their nuclei (see, e.g. 
IComastri et~al]|2002f ). 

3.4. Optical to near-infrared colors 

Figure 10 shows the R— K vs. K-band magnitude (Vega 
magnitudes) for the subsample of objects with ACS 
morphological information (excluding spectroscopically 
confirmed stars). In addition, sources spectroscopically 
identified with BL AGN are marked with green symbols, 
while sources identified with NOT BL AGN are marked 
by yellow symbols. A significant difference in the R— K 
distributions for point-like and extended sources is 
present : while the widths of the two distributions are 
similar [a ^0.8 for both of them), extended sources are 
significantly redder ((R— if }=4.05±0.05) than point-like 
sources ((R — A')=2.91±0.06). When the spectroscopic 
information is also considered, objects with red R — K 
colors are preferentially associated with NOT BL AGN 
(yellow circles), while blue objects arc preferentially 
associated with BL AGN (green circles). 

We have then investigated the distribution of the R — K 
colors as a function of the X-ray hardness ratio (HR), 
a widely used tool to study the general X-ray spectral 
properties of X-ray sources when the nu mber of counts 
is ina dequate to perform a spectral fit. Mainieri et al. 
(|2007t ) (paper IV of this series) have shown that 99% 
of the sources with HR> —0.3 in their subsample can 
be fit by a T= 1.8-2 power-law continuum absorbed by a 
column density (N#) larger than 10 22 cm~ 2 . Conversely, 
sources detected only in the soft band (i.e. HR=-1) are 
most likely unobscured. 

The hardest sources (HR> —0.3, solid histogram in 
the right part of Fig. 10) are mostly associated with 
red and very red objects (R-K>4). This indicates an 
excellent consistency between optical obscuration of the 
nucleus as inferred from ACS (extended morphology), 
optical to near infrared colors, and the presence of X-ray 
obscuration as inferred from the hardness ratio (see also 
Alexander et all 120011; iGiacconi et al.ll2002t IBrusa et all 
20051 [Mainieri et alJl2007t ). Conversely, sources detected 
only in the soft band (i.e. HR=— 1, dashed histogram) 
have preferentially blue R-K colors (> 60% of the 
sources have R — K < 4), ty pical of those of opticall y 
selected, unobscured quasars (|Barkhouse fc HalJl200l . 
It is interesting to note, however, that the observed 
correspondence between hard (soft) X-ray colors and red 
(blue) optical to near infra red colors is not a one-to-one 
correlation (see a lso, e.g.. iGeorgantopoulos et alll2004l ; 
IBrusa et al.ll2005h . In fact, as shown for example by the 
tail toward high R-K colors of the dashed histogram 
in Fig. 10, a non-negligible fraction of red objects is 
associated with very soft X-ray sources. 

3.5. Outliers: mismatch between morphological and 
spectroscopic classification 

As discussed in Sect. 3.1, there is a relatively good agree- 
ment (of the order of 80%) between the morphologi- 
cal and spectroscopic classifications for the subsample of 
sources for which these informations are available. How- 
ever, there are 22 (~ 20%) BL AGN classified as "ex- 
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Absolute B magnitude 

Fig. 11. — Distribution of the absolute B-band magnitude 
for BL AGN (solid histogram, upper panel) and NOT BL AGN 
(solid histogram, lower panel) in the subsample of objects with ACS 
information and spectroscopic classification. In both panels the 
shaded histograms give the contribution of outliers in each class, i.e. 
"extended" BL AGN and "point-like" NOT BL AGN, respectively. 
The dashed line at Mb =-23 marks the "classic" separation between 
Seyfert and Quasars. 

tended" in ACS and 16 (~ 16%) NOT BL AGN classified 
as point-like by ACS. 

The absolute magnitude is an obvious quantity to exam- 
ine in order to try to understand why for these sources 
the morphological and spectroscopic classifications do 
not agree with each other. Figure 11 shows the distribu- 
tion of rest-frame, absolute B-band magnitude for the 
209 extragalactic sources which have both morphologi- 
cal (ACS) and spectroscopic classification, separately for 
BL AGN (upper panel) and NOT BL AGN (lower panel). 
The black filled histograms show the observed distribu- 
tion for the "outliers" , i.e. "extended" BL AGN and 
"point-like" NOT BL AGN, respectively. In order to 
minimize the uncertainties in the K-corrcction, the abso- 
lute B-band magnitudes have been computed from the 
apparent magnitude in the optical filter closest to the 
rest-frame B band at any given redshift. 
All the "extended" BL AGN lie at redshift <1.5, while 
the redshift distribution of the full BL AGN population 
peaks at z~1.5 (sec Figg. 4 and 5). The fact that ACS 
reveals an extended component is likely to be due to 
the fact that the nuclear emission in these objects is not 
dominant with respect to the host galaxy emission, at 
least in the observed ACS band. Indeed, the average 
Mb of the extended BL AGN sources lies in the low-tail 
of the Mb distribution of the entire BL AGN population 
and typically in the Seyfert regime, taking Mb = —23 
as the "classical" separa tion between Seyferts and QSO 
(|Schmidt fc Greed [19831) . 

The unresolved nature of "point-like" NOT BL AGN is 
more puzzling. It can be due, at least in part, to the 
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Fig. 12. — Absolute rest frame B-band magnitude vs. the loga- 
rithm of the full (0.5-10 keV) X— ray luminosity (erg s — 1 ). Symbols 
are as follows: Blue circles = spectroscopically confirmed BL AGN 
classified as point-like from ACS; Yellow circles = spectroscopically 
confirmed BL AGN classified as extended from ACS; Cyan triangles 
= spectroscopically confirmed NOT BL AGN classified as point- 
like from ACS; Red triangles = spectroscopically confirmed NOT 
BL AGN classified as extended from ACS. The horizontal solid line 
marks the "classic" optical transition between Seyfert and Quasars. 



definition we adopted for the segregation of stellar ob- 
jects: especially at high redshift, the ratio between the 
peak flux and the total magnitude cannot be such to 
define them "extended". However, we note that most 
of them occupy the high-luminosity tail of the distribu- 
tion of M B of the overall population of NOT BL AGN, 
and therefore can be classified as high-luminosity, type 2 
quasars from the ava ilable optical spectroscopy (see also 
iMainieri et al.l 120071 for a more detailed analysis of the 
rest-frame X-ray properties). Moreover, they tend to 
be "bluer" (i.e. less absorbed) in the optical to near in- 
frared colors (yellow symbols on blue points in Fig. 10), 
suggesting that we are seeing more directly the emitting 
nucleus. 

Fig. 12 shows the absolute B magnitude as a function of 
the .5-10 keV luminosity (as computed in lMainieri et al.l 
[2007j) for the subsample of sources with spectroscopic 
redshift and good X-ray counting statistics. Circles 
mark BL AGN (blue for the point-like and yellow for 
the extended sources) while triangles mark NOT BL 
AGN (cyan for the point-like and red for the extended 
sources). Sources with Lx < 10 42 erg s" 1 are mostly 
associated with bright (I<18), extended objects (red 
triangles) at z<0.2; two of them are detected only in 
the soft band. Among the two sources detected also at 
higher energies, there is source the can d idate Compton 
Thick AGN discussed in lHasinger et all (j2007T l. 
The region at luminosities higher than 10 44 erg s" 1 
is mainly populated by point-like, BL AGN (blue 
circles). Only 8 sources classified as NOT BL AGN lie 



in this part of the diagram: the high X- ray luminosity 
classi fies them as candidate Type 2 QSO ( Maimeri et alJ 
l2007t ). Conversely, the majority of the sources in the 
middle part of the X-ray luminosity region are associ- 
ated with extended, NOT BL AGN (red triangles) and 
low redshifts BL AGN (Seyfert 1 galaxies, yellow circles). 

4. DISCUSSION 

Comparing the optical color-color properties of AGNs 
in the COSMOS field with those of field objects (see 
Fig. 8), we estimate that X-ray data with a flux limit 
of So.5_2keV ~ 10~ 15 erg cm -2 s -1 recover at least half 
(50%) of the AGN candidates which would be selected 
on the bas is of their ultraviolet excess (U-B<0.3) at 
B<23 (see IZamorani et alJ 119991 for a similar estimate 
at a somewhat brighter X-ray flux from R OS A T obser- 
vations in the Marano field). This fraction of 50% has 
been obtained assuming that all objects with ultraviolet 
excess are AGN. Since at these magnitudes blue stars 
and compact galaxies may represent ~ 30% of the 
samples of point-like objec ts with ultraviolet excess 
(jMignoli fc Zamorani I l2001h . we conclude that the 
efficiency of present X-ray data in detecting AGN with 
ultraviolet excess at B<23 is likely to be as high as 70%. 
We note that the large majority of the BL AGN selected 
via the UV excess method have MB< —23 (see Fig. 11) 
and are therefore cl assified as quasars, con firming previ- 
ous results (see e.g IZamorani et al.1 Il999l and reference 
therein). We also note, however, that the photometric 
catalog we use has deeper limiting magnitudes with 
respect to all the previous studies on optically selected 
UV-excess QSO and therefore we start to explore and 
pick up lower-luminosity (e.g. Seyfert-like) objects up to 
z~ 1.5. On the other hand, a significant fraction (~40%) 
of the X-ray selected AGNs, especially those without 
broad lines in their optical spectra, would have not been 
easily selected as AGN candidates on the basis of purely 
optical criteria, either because of colors similar to those 
of normal stars or field galaxies at z~l-3, or because of 
morphological classification not consistent with that of 
point-like sources (see Fig. 10). These results highlight 
the need for a full multiwavelength coverage to properly 
study and characterize the AGN population as a whole. 
With the current spectroscopic coverage (~38%) the 
observed redshift distribution and spectroscopic clas- 
sification of identified X-ray sources could be biased 
against low-redshi ft BL AGN and hig h-redshift NL AGN 
(see Sect. 2.5 and lTrump et al.ll2007l for a more detailed 
disc ussion on AGN sele ction effects in COSMOS; see 
also iTreister et al.1 l2004j ). However, according to the 
most recent mod eling of the X-ray luminosity function 
evolution (e.g. lUeda et all l2003t iBarger et all 120051; 



Hasinger. Mivaii fc Schmidt! l2005t lLa Franca et al.l 
20051 ). the paucity of low-redshift BL AGN and high- 



rcdshift obscured NL AGN may be at least partly 
real and due to a luminosity dependent evolutionary 
behaviour of X-ray selected sources. The space density 
of high luminosity (Lx > 10 44 erg s _1 ) quasars decreases 
steeply below z~l, while the decrease of the space 
density of lower luminosity objects in the s ame redshift 
range is much slower (see, e.g., Fig. 9 in iFiore et aTl 
|2003|) . Moreover, there are now increasing evidences 
that the relative ratio between obscured and unobscurcd 
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AGN is a d e creasi ng function of the X-ray luminosity 
(lUeda et all 120031: lHasingerl 120041 : lLa Franca et all 
120051: lAkvlal et al.1 |2006| ). As a consequence, at low 
rcdshifts the X-ray source population is dominated by 
low luminosity, often obscured AGN, while at higher 
rcdshifts luminous unobscured quasars take over. This is 
in qualitative agreement with the observed change of the 
Type 2/ Type 1 ratio (R) in the present sample when 
the typical Seyferts (32/33, R-l at L x ~ 10 42 - 10 44 
erg s _1 ) and Quasars (8/44, R-0.18 at L x > 10 44 erg 
s _1 ) luminosities are considered (see Fig. 12). 

At a limiting flux of 2xl0 -15 erg cm~ 2 s _1 in the soft 
band, where the sky coverage has decreased to 50% of 
the total 1.3 deg 2 area covered by the 12F pointings, no 
object at z>3 is present so far in our current spectro- 
scopic sample. On the one hand, the lack of high redshift 
sources could be due, at least in part, to the still limited 
spectroscopic completeness, which is of the order of 50% 
at Iab<22 and about 25% at fainter magnitudes, with 
no spectrosco pi c reds hifts for Iab>24. As a comparison, 
iMurrav et all lj2005ft have 14 z> 3 quasars in their 
9 deg 2 X-Bootes survey, at a limiting flux about one 
order of magnitude brighter than our survey and with 
a similar spectroscopic completeness. By rescaling the 
area and the limiting fluxes, this would translate in ~ 3 
z> 3 quasars in the XMM- COSMOS 12F area, that 
is still consistent with zero being observed. On the 
other hand, the number of AGN at z>3 expected from 
popu l ation synthesis models (see also Hasinge r et alJ 
I2007t iGilli. Comastri fc Hasingeil l2007t ) is ~30 and 
a random sampling at ~25% completeness, to take 
into account the current follow-up spectroscopy, would 
therefore predict about 7 quasars at z>3, which are not 
observed. The forthcoming spectroscopic observations 
will increase the spectroscopic completeness, allowing to 
verify if a sizable number of high redshift (z>3) quasars 
is present among the optically faintest counterparts 
and among the optically unidentified X-ray sources, 
or if some of the model assumptions, mainly based on 
extrapolations of the lower redshift (z< 3) luminosity 
functions, should be refined. 

5. SUMMARY 

In this paper we presented the optical identifications for 
a large subsample (~700) of X-ray sources detected in 
the first 1.3 deg 2 of the XMM-COSMOS observations, 
down to a 0.5-2 keV limiting flux of ~ 10~ 15 erg cm~ 2 
s . The X-ray counterparts have been identified in 
optical (I-band) and near infrared (K-band) catalogs, 
making use of the "maximum likelihood" technique. 
The combined use of the two different catalogs allowed 
us to test the identification procedure and turned 
out to be extremely useful both for isolating input 
catalogs problems and for identifying optically faint 
counterparts. Overall, 90% (626) of the X-ray sources 
have been unambiguously associated with optical or 
near-infrared counterparts. Twenty eight X-ray sources 
have 2 possible optical counterparts with comparable 
likelihood of being the correct identification. These 
sources have been classified as "ambiguous" and we 
tentatively identify them with the optical counterpart 
with the highest LR value. The sample of proposed 



identifications therefore comprises 654 objects, while 
for the remaining 41 X-ray sources it was not possible 
to assign a candidate counterpart. We will use the 
forthcoming Chandra data (granted in A08) to defini- 
tively discriminate the ambiguous and false associations 
and we also predict that a large fraction of these very 
faint, unidentified objects will be resolved with the 
inclusion of Spitzer (IRAC and MIPS) data in the source 
identification. 

We then cross-correlated our proposed opt ical counter- 
parts with the Subaru multi-co lor catalog (ICapak et al.l 
l2007h . the HST/ACS data (jLeauthaud et all l2007h 
and the first results from the m assive spectroscopic 
campaigns in th e COSMOS field (jTrump et all 120071 ; 
iLillv et alll2007l ). Our analysis reveals that for ~80% 
of the X-ray sources counterparts with spectroscopic 
rcdshifts (a total of 248) there is a good agreement be- 
tween the spectroscopic classification, the morphological 
parameters, and optical to near infrared colors: the large 
majority of spectroscopically identified BL AGN have a 
point-like morphology on ACS data (see Sect. 3.1), blue 
optical colors in color-color diagrams (see Sect. 3.2), and 
an X-ray to optical flux ratio typical of optically selected 
quasars (see Sect. 3.3). Conversely, NOT BL AGN are 
on average associated with extended optical sources, 
have significantly redder optical to near infrared colors 
and span a larger range of X-ray to optical flux ratios 
(see Sect. 3.3 and 3.4). When the X-ray information is 
also considered, we found that hard X-ray sources are 
preferentially associated with extended s ources and are 
"redd ish" (~60% with R-K>4, see also iMainieri et~al"1 
l200l . while sources detected only in the soft band 
are mostly associated with point-like objects and have 
"blue" optical to near-infrared colors. 
Comparing the optical color-color properties of AGNs 
in the COSMOS field with those of field objects (see 
Fig. 8), we estimate that X-ray data with a flux limit 
of So.5-2kcV ~ 10 -15 erg cm" 2 s" 1 recover ~ 70% of 
the AGN selected on the basis of their ultraviolet excess 
at B<23. This fraction will rise up to (90-100)% at 
So 5-2koV ~ 5 x 10~ 16 erg cm~ 2 s -1 , the final depth of 
the XMM-COSMOS survey. 

About 20% of the sources show an apparent mismatch 
between the morphological and spectroscopic classifi- 
cations. Our analysis indicates that, at least for BL 
AGN, the observed differences are largely explained by 
the location of these objects in the redshift-luminosity 
plane (see sect. 3.5). Our analysis also suggests a change 
of the Type 2/Type 1 ratio as a function of the X-ray 
luminosity, in qualitative agree ment with the results 
from X-ray spectral analysis (|Mainieri et al.l l2007f ) 
and the most recent modeling of the X-ray luminosity 
function evolution (Figg. 5 and 12). 
Although the Magellan/IMACS and VIMOS/zCOSMOS 
spectroscopic campaigns will continue to obtain rcdshifts 
for AGN in the COSMOS field, we expect that a fraction 
of the order of 30% of the sources will not have spec- 
troscopic redshifts, due to the faintness of their optical 
counterparts (~ 15%) and the limi tation of multi-sli t 
spectroscopy due to mask efficiency (jlmpev et al.ll2007t ) . 
We have tested a photometric code specifically designed 
to obtain redshifts for X-ray selected sources on the sub- 
sample of sources with spectroscopic data, and showed 
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that it is possible to obtain a reasonable estimate of the 
sources redshifts (\z spec - z p hot\ < 0.15 x (1 + z spC c)), for 
~70% of the sources. We plan to use the same code to 
derive the redshifts for the faintest X-ray counterparts. 

Summarizing, we were able to perform a combined mul- 
ticolor, spectroscopic and morphological analysis on a 
statistical and meaningful sample of X-ray sources, and 
this work represents one of the most comprehensive mul- 
tiwavclcngth studies to date of the sources resp onsible 
of most of the XRB (see also lEckart et al.ll2006D . This 
work is just the first phase of COSMOS AGN studies; 
the results presented in this paper clearly show the need 
for a full multiwavelength coverage to properly study and 
characterize the AGN population as a whole. As more 
data will become available, the selection of COSMOS 
AGNs by all available means - X-ray, UV, optical, near- 
IR, and radio - will build up to give the first bolomet- 
ric selected AGN sample, fulfilling the promise of many 
years of multi-wavelength studies of quasars. 
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